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A PROGRAM  TO  PLOT  A TRACK  AND  BATHYMETRY  OR  MAGNETIC 
PROFILE  ON  A POLAR  STEREOGRAPHIC  PROJECTION 


1.0  IDENTIFICATION 

1.1  Title 

Program  to  Plot  a Track  and  Bathymetry  or  Magnetic  Profile  on  a Polar 
Stereographic  Projection. 

1 .2  Identification  Name 
Track. 

1.3  Classification  Code 
None. 

1.4  RCC  Identification  Number 
None. 

1.5  Entry  Points 
TRACK. 

1.6  Programming  Language 
Language:  CDC  3600/3800  Fortran. 

Routine  Type:  Program. 

Operating  System:  Drum  Scope  2.1. 

1.7  Computer  and  Configuration 
CDC  3800. 

1.8  Contributor  or  Programmer 

Marilyn  L.  Blodgett,  Code  8176MB,  Long  Range  Propagation  Section, 
written  for  the  Environmental  Sciences  Section,  Acoustics  Division. 

.1 .9  Contributing  Organization 

NRL  - Naval  Research  Laboratory,  Washington,  D.C.  20375. 

1.10  Program  Availability 

If  supplied  with  a magnetic  tape,  the  Environmental  Sciences  Section. 
Acoustics  Division,  will  make  a copy  of  this  program. 

1.11  Verification 

This  program  has  been  used  and  tested  by  the  Environmental  Sciences 
Section,  Acoustics  Division,  for  several  months. 

1.12  Date 
April  1976 


Manuscript  submitted  September  3,  1075. 
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PURPOSE 

✓ 

2.1  Description  of  the  Routine 

This  program  reads  the  data  collected  by  an  oceanographic  or  geophysical 
experiment  from  a magnetic  tape  and  plots  the  track  and  bathymetric  or 
magnetic  value  perpendicular  to  the  track  as  a profile.  We  use  the  format 
recommended  by  the  National  Research  Council  of  the  National  Academy 
of  Sciences  with  one  slight  modification  for  the  input  data  tape.  There  is 
one  logical  record  (of  80  characters)  for  each  data  point.  The  different 
types  of  data  (bathymetry  and  magnetics)  are  separated  by  an  end-of-file 
mark  with  a double  end-of-file  mark  at  the  end  of  all  the  data. 

Before  the  program  reads  this  input  tape,  it  reads  two  cards.  The  first  card 
defines  the  actual  data  format  on  the  input  tape  (the  format  varies  for  the 
two  types  of  data).  The  second  card  specifies  the  number  of  files  to  be 
skipper,  over  on  the  first  input  tape,  the  physical  height  of  the  map  to  be 
drav'.i,  the  actual  latitude  and  longitude  values  to  be  included  on  the  grid, 
tlv-  dates  of  the  data  on  the  first  input  tape  to  be  considered  for  plotting, 
the  actual  values  to  be  plotted,  and  the  units  per  inch  for  plotting  the 
bathymetric  or  magnetic  profiles  along  the  track. 

With  all  the  required  parameters  defined,  the  program  starts  to  read  the  in- 
put tape  one  record  at  a time.  Each  record  is  checked  to  see  that  the  fix 
falls  on  the  defined  grid  and  that  it  was  taken  on  or  between  the  two  spec- 
ified dates.  Only  those  points  which  meet  both  requirements  are  stored  in 
core.  The  program  continues  reading  the  first  input  tape  until  it  reads  an 
end-of-file  mark  or  a fix  taken  after  the  last  specified  date.  If  there  are 
additional  input  tapes,  the  program  reads  them  in  a similar  manner.  The 
beginning  and  end  dates  for  each  new  input  tape  are  contained  on  an  Extra 
card.  A maximum  of  four  input  tapes  can  be  used.  When  all  the  input 
tapes  have  been  read,  the  program  prepares  to  plot  the  track  and  the  speci- 
fied values,  either  bathymetry  or  magnetics. 


The  track  is  plotted  on  a polar  stereographic  projection  which  is  drawn 
exactly  to  scale.  The  grid  may  be  blown  up  to  any  reasonable  size.  The 
largest  grid  we  have  defined  is  1 degree  of  latitude  equals  20  inches.  The 
number  of  degrees  of  longitude  included  in  the  grid  will  depend  on  the 
scale  of  the  entire  grid  and  the  specific  area  of  interest.  In  the  case  of  1 
degree  of  latitude  equals  20  inches,  no  more  than  10  degrees  of  longitude 
can  be  included  in  the  grid.  Since  the  projection  is  drawn  exactly  to  scale, 
a mosaic  can  later  be  built  of  the  entire  area.  Depending  on  the  type  of 
data  read,  the  profiling  values  will  be  either  uncorrected  fathoms,  uncor- 
rected meters,  or  residual  magnetic  intensity. 

2.1.1  Bathymetry  Data 

The  program  reads  the  year,  date  (month  and  day),  hour,  minute,  latitude, 
longitude,  and  uncorrected  fathoms  from  the  input  tape  according  to  the 
specified  format.  The  southern  latitudes  and  the  western  longitudes  are 
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preceded  by  a negative  sign.  The  program  can  convert  uncorrected  fathoms 
to  uncorrected  meters.  The  track  is  plotted  in  a continuous  straight  line, 
and  the  profiling  series  is  either  uncorrected  fathoms  or  meters  multiplied 
by  -1  to  drop  it  below  the  track. 

2.1.2  Magnetic  Data 

The  program  reads  the  year,  date  (month  and  day),  hour,  minute,  latitude, 
longitude,  and  residual  magnetic  intensity  from  the  input  tape  according 
to  the  specified  format.  The  southern  latitudes  and  the  western  longitudes 
are  preceded  by  a negative  sign.  The  track  is  plotted  in  a continuous 
straight  line,  and  the  profiling  series  is  residual  magnetic  intensity. 

2.2  Problem  Background 

Program  Track  was  written  so  that  the  researcher  can  build  a profile,  either 
magnetic  or  bathymetric,  along  the  track  from  which  the  data  were  taken. 
Presenting  data  in  this  manner  will  show  bathymetric  or  magnetic  trends  in 
relation  to  the  geographic  area. 


USAGE 

3.1  Calling  Sequence  or  Operation  Procedure 
Not  applicable. 

3.2  Arguments,  Parameters,  and/or  Initial  Conditions 
Not  applicable. 

3.3  Space  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage 

5127  octal  (2647  decimal)  locations  exclusive  of  system  library 
functions. 

3.3.2  Common  Blocks 
Blank  common 

/l/,  131,  1 5/,  HI,  181,  19/,  1 10/. 

3.3.3  Temporary  Storage 
None. 

3.4  Messages  and  Instructions  to  the  Operator 
None. 

3.5  Error  Return,  Messages,  and  Codes 
None. 

3 .6  Informative  Messages  to  the  User 
None. 
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Input 

Tho  actual  format  of  the  data  on  the  input  tape,  the  map  parameters,  and 
the  command  words  are  read  in  via  input  cards.  The  track  and  the  data 
to  be  profiled  are  read  in  via  magnetic  tape  on  logical  units  15  through  18. 
Apendix  A presents  samples  of  our  data  formats  on  the  input  tape.  Appen- 
dix B is  a complete  description  of  the  input  deck  setup. 

Output 

The  program  prints  on  the  standard  piinter  (logical  unit  61)  the  data  format, 
chart  parameters,  number  of  data  points  read  in,  and  the  number  of  data 
points  plotted  on  the  map  for  both  the  track  and  the  profiled  data.  Appen- 
dix C presents  sample  profiles,  and  Appendix  D presents  a sample  output 
listing.  The  program  writes  the  plotting  instructions  on  a magnetic  tape 
(logical  unit  40). 

Formats 

Appendix  B describes  the  program  deck  structure. 

External  Routines  and  Symbols 

ATAN2,  SQRTF,  SINF,  COSF,  ATANF,  SPACE00,  BACKFILE  SKIPFILE 
PLOTS,  NUMBER,  STOPPLOT,  PLOT,  SYMBOL. 

Timing 

The  time  required  depends  on  the  size  of  the  grid  and  the  number  of  data 
read  and  plotted. 

Accuracy 

The  grid  is  reproduced  exactly  to  scale. 

Cautions  to  Users 
None. 

Program  Deck  Structure 

Appendix  B describes  the  program  deck  structure. 

References  - Literature 

R.L.  Parker,  “The  UCSD  Hypermap  Programs,”  University  of  California 
Sar  Diego. 

M.J.  Kertyzak  and  J.D.  Phillips,  “GRENFY,”  Woods  Hole  Oceanographic 
Institute,  Woods  Hole,  Massachusetts. 

M L-  Blodgett  and  J.V.  Massingill,  “A  Program  for  Storing  Oceanographic 
Data  on  Magnetic  Tape,”  NRL  Report  7861,  March  1975. 
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METHOD  OR  ALGORITHM 
Not  Applicable. 

FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 
The  flow  chart  and  listing  are  given  in  Appendixes  E and  F. 

COMPARISON 

No  other  known  programs  are  available  for  comparison. 

TEST  METHOD  AND  RESULTS 

The  program  has  been  used  and  tested  successfully  on  a Calcoinp  plotter. 

REMARKS 

None. 
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APPENDIX  B 

Deck  Assembly  for  Program  Track 


21  End  of  File  card 
20  Stop  card 
19  Continue  card4 

18  Plot  Anoml  card 
17  Special  Anoml  card 

16  Plot  Value  Card 

15  Extra  card  (one  for  each  additional  input 

tape) 

14  Special  Value  card 
13  Actual  Title  card3 
12  Title  card3 


^The  program  uses  scratch  tapes  on  logical  units  20  and  05,  but  no  Equip  cards  are  required,  since  the 
drum  is  used. 

^If  the  Fortran  source  deck  is  used  instead  of  the  binary  deck,  a Fortran  card  is  required  after  the  Bank 
card.  In  addition,  a Scope  card  and  Load  card  must  follow  the  source  deck. 

3Thcse  two  cards  are  not  required  by  the  program;  both  cards  may  be  present  or  both  omitted. 

4This  card  is  used  only  if  another  plot  is  desired.  It  is  to  be  followed  by  a second  set  of  input  cards 
(1018). 
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Card  litle 
Job 


Equip 


Bank 


Program 


Run 


Actual 

Format 


Column 

Number 

1-21 


1-13 


Deck  of 
Track 


1-13 


1-? 


Description 

7/9  JOB,  Charge  No..  ID  No.,  time. 
See  page  2-2  of  the  3600/3800  Com- 
puter System  Drum  Scope  Manual. 

7/9  EQUIP,  40=**,  WO,  LO 
7/9  EQUIP,  15=**,  RO,  HI 
7/9  EQUIP,  16=**,  RO,  HI 
7/9  EQUIP,  17=**,  RO,  HI 
7/9  EQUIP,  18=**,  RO,  HI 
40,  15,  16,  17,  18  = logical  unit 
numbers. 

RO  = read  only.  WO  = write  only. 
LO  = low  density.  HI  = high  density. 

-/0/7/9  BANK,  (0),  /l/  See  page  7-17 
of  the  3600/3800  Computer  System 
Drum  Scope  Manual. 

This  is  the  main  program  with  asso- 
ciated subroutines.  If  the  Fortran 
source  deck  is  used  instead  of  the 
binary  deck,  a Fortran  card  is  re- 
quired after  the  Bank  card.  The 
Fortran  card  reads  7/9  FTN,  L,  R,  X. 
In  addition  a Scope  card  with  SCOPE 
starting  in  column  10  and  a Load 
card  must  follow  the  source  deck. 

7/9  RUN,  T,  P,  R,  M,  D 
T = time  limit  in  minutes. 

P = Maximum  number  of  print  or 
write  operations. 

R,  M,  D may  be  left  blank.  See  page 
2-15  of  the  3600/3800  Computer 
System  Drum  Scope  Manual. 

(13X,  12,  14,  1X12,  F3.1,  F8.4, 

F9.4,  28XF5,  5X)  This  format  should 
be  replaced  by  the  desired  input  for- 
mat. The  format  must  be  enclosed  in 
parentheses  and  left-justified.  Via 
this  format  the  program  reads  the 
year,  date,  hour,  minute,  latitude, 
longitude,  and  value  for  the  profiling 
series  (uncorrected  fathoms  for  ba- 
thymetry and  residual  magnetic  inten- 
sity for  magnetics). 
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Number 


11 


Column 

Card  Title  Number 


Map  Parameter  2 


4 


6 


9-10 


11-12 


Description 

The  formats  for  reading  the  two 
data  types  on  our  input  tapes  are: 

• Bathymetry 
(13X12,14,1X12,  F3.1, 

F8.4,  F9.4,  10XF5.1) 

• Magnetics 

(13X12, 14, 1X12,  F3.1, 

F8.4,  F9.4,  28XF5). 

-1,  0,  or  1 

-1  = multiply  uncorrected  fathoms 
by  -1  to  drop  the  value  series  below 
the  track. 

0 = plot  the  anomaly  value  as  read 
from  the  input  tape.  This  parameter 
is  used  to  plot  the  residual  magnetic 
intensity. 

1 = convert  uncorrected  fathoms  to 
uncorrected  meters  and  multiply  by 
-1  to  drop  the  profiling  series  below 
the  track. 

0 or  1 

0 = plot  only  the  track. 

1 = plot  both  the  track  and  the 
profiling  series. 

0 or  1 

0 = plot  all  data  which  falls  on  the 
defined  grid. 

1 = plot  all  data  which  falls  between 
the  southern  degree  of  latitude  plus 
one  degree  and  the  northern  degree 
of  latitude. 

2 

Number  of  degrees  between  latitude 
lines  drawn  on  the  grid. 

1 

Number  of  degrees  between  the 
longitude  lines  drawn  on  the  grid. 
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Number 


Card  Title 


Column 

Number 

13-14 


15-22 


23-30 


31-38 


39-46 


47-54 


55-62 

64 


Description 

0,  1,  or  2 

Number  of  files  to  be  skipped  over 
our  input  tape.  There  are  a maxi- 
mum of  three  files  on  our  Geodata 
tapes. 

1000.00 

Units  per  inch  for  plotting  the  anom- 
aly along  the  track.  The  maximum 
is  2000  gammas  per  inch  on  the  map 
surface.  A value  of  1000  means 
that  a profiling  value  of  1000  gam- 
mas would  be  plotted  1 inch  above 
the  track.  The  remainder  of  the 
anomaly  data  would  be  scaled  ac- 
cordingly. 

20.0 

Physical  height  of  the  chart  to  be 
drawn.  To  obtain  this  figure,  you 
must  measure  the  actual  physical 
height  from  an  existing  map. 

82.5 

The  degree  of  latitude  at  the  base 
of  the  chart  (the  southernmost  lat- 
itude). This  value  may  be  either  a 
whole  or  a half  degree.  (Southern 
latitudes  are  preceded  by  a minus 
sign.) 

84.5 

The  northernmost  degree  of  latitude. 
The  difference  between  the  degrees 
of  latitude  should  be  an  integer. 

-15.0 

The  westernmost  degree  of  longitude. 
(Western  longitudes  are  preceded  by 
a minus  sign.) 

05.0 

The  easternmost  degree  of  longitude. 

1,  2,  3,  or  4 

Number  of  input  tapes,  with  the 
maximum  being  four  tapes. 
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Column 

Number  Card  Title  Number 


65-72 


73-80 


12  Title 


1-5 


13  Actual  Title  1-80 


14  Special  Values  1-14 


15  Extra  1-4 


Description 


02251600 

The  date  and  time  of  the  first  data 
point  to  be  plotted  from  the  first 
input  tape.  Columns  65-66  = month, 
67-68  = day,  69-70  = hours,  and  71- 
72  = minutes. 

02280830 

The  date  and  time  of  the  last  data 
point  to  be  plotted  from  the  first 
input  tape.  All  data  taken  on  and 
between  the  dates  and  times  of  the 
first  and  last  data  points  will  be 
plotted  if  they  fall  within  the  de- 
fined chart. 

TITLE 

This  command  allows  the  user  to 
label  the  chart.  This  is  a non- 
obligatory  card. 

ARCTIC  BASIN 

The  appropriate  title  may  be  punched 
anywhere  in  the  80  columns.  This 
is  a nonobligatory  card. 

SPECIAL  VALUES 
This  command  allows  the  user  to 
associate  the  name.  Values  with  the 
series  of  data  points  read  from  the 
input  tape(s).  The  program  will 
store  only  those  data  points  which 
fall  on  the  defined  chart  and  which 
were  taken  on  or  between  the  two 
dates  specified. 

0,  1,  or  2 

Number  of  files  to  be  skipped  over 
on  the  second  input  tape.  There 
must  be  an  Extra  card  for  each  ad- 
ditional input  tape.  Since  there  is 
a maximum  of  four  input  tapes,  the 
maximum  number  of  Extra  cards  is 
three. 
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Number 


16 


17 


18 


19 


20 


21 


Column 

Card  Title  Number 

5-12 


13-20 


Plot  Values  1-11 


Special  Anoml  1-13 


Plot  Anoml  1-10 


Continue 


Stop 


End  of  File 


Description 


02251600 

The  date  and  time  of  the  first  data 
point  to  be  plotted  from  the  second 
input  tape.  The  dates  for  the  first 
input  tape  are  on  the  Map  Param- 
eter card. 

02280830 

Date  and  time  of  the  last  data  point 
to  be  read  and  plotted  from  the 
second  input  tape. 

PLOT  VALUES 

This  command  causes  the  named 
series  to  be  plotted. 

SPECIAL  ANOM1 
This  command  allows  the  user  to 
plot  the  profiling  series.  Use  only  if 
there  is  a 1 in  Column  4 of  Card 
No.  11. 

PLOT  ANOM1 

This  command  causes  the  profiling 
series  to  be  plotted.  Use  only  if 
there  is  a 1 in  Column  4 of  Card 
No.  11. 

This  card  is  used  only  if  another 
plot  is  desired.  It  should  be  fol- 
lowed by  a set  of  control  cards 
(cards  10  through  18).  The  pro- 
gram will  not  rewind  the  input  tapes. 
It  will  continue  reading  where  it 
left  off  unless  told  to  skip  to  another 
file  by  the  Map  Parameter  card. 

STOP 

This  command  terminates  the 
program. 

Terminates  the  run. 
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DATA  FORMAT  . . . (13X , 12 , 14 , 1X12 

CHART  PARAMETERS 

SOUTHMOST  LATITUDE  74.0 

WESTMOST  LONGITUDE  0.0 


PROGRAM 

READ  IN 

1279 

POINTS 

PROGRAM 

PLOTTED 

1278 

POINTS 

PROGRAM 

READ  IN 

1279 

POINTS 

PROGRAM 

PLOTTED 

975 

POINTS 

,F3.1,F8.4,F9.4, 10XF5 . 1 , 23X) 

NORTHMOST  LATITUDE 
EASTMOST  LONGITUDE 

ON  THE  MAP 
ON  THE  MAP 


76.0 

10.0 
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APPENDIX  E 
Flow  Chart 


APPENDIX  F 
Source  Language  Listing 


C 

C 


PRGqRam  Track 
DIMENSION  I0LH  254  ) 

DIMENSION  Name<6)  , IPM20) 

real  latmjn.latmax.conomin.longpax  , la t 
COMMON  DEL*!. tEL®N.»FElAT,XPELPN 
C£MMGN  LX 

C6MM0N  H.DIST.ANOPCK.CFANGE.NP.LKK.KNUM, I NUM » G INCH 
COMMON  ANOM(  2) . I 1 I »J„J,KKK,XLAST 
COMMON  LAY ( 2 ) 

COMMON  inp, jo.u.v.xlat.xlon 

COMMON  P0L at (POLONC . POT  ilMlN.LPAX.VMlN, VPAX,HEIGHT,NLAT,NL0N, 
♦IPR0J, IB0X, SCALE. ISYP8, ILINE 

COMMON  NALL, I ETaRT (li ),NAME5 <10,6 >,leng, last 

COMMON  TES,TITL£(1C),1C0L,IA(7«) 

COMmON/I/JDaTE (2)  , HE  DN  ( 2 ) 

CGMmON/J/LATPIN.LATMAX 
COmmON/5/JUDYI , JUCV2, 1TM1, JTP2 

COMmON/7XLONGP in.longmax 
common/bx i r m 

COMMON/lO/IEXTRA,  1SK1P 


DATa(kEY1*4HC0NT),  (KEY 2*4HREA£),(KEy3*4HC6MP>,(KEY4«4HSPEC),  (KEYS* 
13PMAP,.<KEY6MHT1TL,.<XEY7.4EhLND),,KfcYB.4HPLBT;.«KEY9MMpiilS);;!E 
2Y1O«4hST0P),(KEY11»4PEFAS),(KEY12*4HST0R),( KEY13«4MRECA ) . (KEY14MH 
3STHE)l(KEY15.4HP0LE),(KEYl«.4hTLRN>,(K6Yl7.4HFPLe) 

REWIND  15  fc 
6601  REW|ND  20 


KNUM*0 
INUM*0 
JVUvO 
NALL«0 
LAST*0 
rewind  05 
REWIND  06 
I NP*60 
10*61 

ISTARTU)«1 

YES»-100.0 

IC0L«1OO 

KMK  IS  A INDICATOR  [F«0  WILL 
ir»l  WILL  CALCULATE  ANOMALY 

KKKil 
D I S T ■ 1 0 , 

AN0MCK=15OO, 

CMnGE*20, 


8NLT  READ  TRACK  SERIES  FR8N  CAICM 
SERIES  AS  WELL 


CALL  PLOTS ( IPLP ,254.40(29) 

Rt*D  ( 1 NP, 1000  ) IFP 
WRITE* 10.1001 ) JFP 
1000  FORMAT ( 2 o A 4 ) 

1001  F0RMaT(///*  DATA  FORMAT  tj.*,20A4) 

SCAlE.O, 

W I DTH«Hfc I GHT 


Dl  701  |Xali |SK|P 


19 


BLODGETT  AND  MASSJNGILL 


701  CALI,  SKIPHLEUS) 

CIST-10 ,0 
CH*NGfc*20,0 
AN0MCK  *15Q0 , 0 
lUNEfl 
700  I PR@J«7 
peigNQ»o  .0 
PeiATi90,0 
F * i 017*933 
SKAlE»C0S(LATP'A|*F  > 

DEL9N«L0NOMAX-t,SNGMA 

delat^latm4x-l/ '■'•IK 

Xl,0N»O  5 
XL*T"LATMAX 
CALL  CSNV(XLAT) 

VP  A X A V 

XLAT^LATMJN 
CALL  C0NV(XL*T) 

VPIUV 

UP. Ax«  (3, 1415926536/36 0,  >*DELCN»SKALE 
UP I N*  »UMAX 
910  WRITE'  1 0.10  03  ) 

1003  F0RMT  ( lHQilfthCHAHT  FARAPETfeRS) 

WRITE? IO#1004)LATPIN,LATPAX 

100«  FCRM*  T ( 1H  ,2OhS0UTPPeST  LATITUDE  ».F10 . 1 . 10X . 2OHN0RTHM0ST  LATITUDE 
1 iFIO  • 1 ) 

WR 1 TE ( 1 0 . 1003  > L0PCPIN.LENOPAX 

1009  FERHATdH  »2OHWESTP05T  LING  1 TIDE  •ri0.1.1OX.2OHEA«TM0ST  LINOITUDE 
1 iFIO • 1 ) 

10  CALL  NEXT(KET,KAME) 

IF  (KEv.EO.KETl)  G£  T { £600 
IF (KEY.EO.KET4)  oe  Te  J055 
IF (KEY.EO.KET6)  oe  TE  600 
IF (KEy.EO.KEY8)  CALL  8LTPUT < N APE , 1 ) 

IF  (KEY  .EQ,  KEY10)  G0  TC  100 
GE  TG  10 

1059  CALL  0)HER(LAT<1>.LAT(LAST*1> J 

LAST«L*ST  ♦ NF 

IF ( LK K .E0.2)  GC  T£  1C 

WR I TE ( 10.1050)  N P 

1090  FERMAT (16H0PRCGRAH  REAC  IN* I 9 » 2X6HP0 J NT$ ) 

GE  TG  10 

600  READ( JNP.8000)  title 
8000  F 6Rm*T ( 1 C A 0 ) 

WRITECI0.6OOO)  TITLE 
6000  FERMAT(1H0,12PMAP  TITLE  A 1PA8  I 
TES»100  i 0 
GE  T0  10 

6600  CALL  PL0T (WIDTH  *10,0#0,0i»3) 

GC  TB  6601 
100  CALL  PL0TS(O,O) 

CALL  6TBPPLBT 
END 


20 


NRL  REPORT  7929 


program  lenqth 

01166 

entry  pe i nts 

track 

0 C 5 4 1 

BLOCK  NAMES 

i 

00315 

00004 

3 

00  002 

5 

00004 

7 

0 0002 

8 

00024 

internal  svmbbls 

10 

00002 

OflOENTRy 

THENO, 

OQCC IC  T , 

plots 

SKIFF IUE 

C0Ny 

NEXT 

output 

OTHER 

PLOT 

STOPPlOT 

COST 


REH, 


00233  SVMBSLJ 


TSN 

STH 


» 

I 


SL8, 

SI!  i 

onsjkgl, 


1 CENT 


TRACK 


C 

ccc 

c 


SUBROUTINE  CGNV  <LAT » IPEN, IFRSv ) 

BASIC  SUBROUTINE  CONTAINING  12  STANDARD 

real  lat 

DIMENSION 
DIMENSION 
DIMENSION 

dimension 
CEMMON 
CEMMON 

cemmon 
common 

C6MMGN 
CEMMGN 

cgmmgn 

1 i i L JNe 

CGMmON 
cgmmgn 


PROJECT  1 0NS i 

ME *C ( 2 ) 

S(2),P(2J 

lat  < 2 ) 

A<4) 

DELAT,CEL0N,XPGLAT,XPeLON 

LX 

lllMni?T'*N(:,f'CK'CMNGE,hF  -U'l'tKNUM,  INUR.OINCH 
AN0M(  2)iIIII»vJv|KKK»XLAST 
X ( 2 ) 

INP,I®,U(2),XLAT,XLGN 

POL AT, FRJ0,FGT,LV<4), HEIGHT, NL AT, NL0N.JI I ( 2 > , SCALE » ISYMB 


<11J  »n*ME5  (10»6)fLENG,LAST 
TES,TITLE»J0MCBL,IA(76I 

CGMmON/1/|DATE<2)  , HEEM2  ) 
C6Mm6n/9/UD|FF,DIFF , III  IN,  VV  JN 
DAT^IRAT 10*1,00092), (ECCSQ*0,0087227) 
data (F», 0174533 >«  CLLAST»9999,  ) 
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RCF»37. 2957795 
I C A Y ■ 0 • 0 
JPKM*o 
2£R0»O.O 
ZN I NE«99 . 0 
PJ"3, 14159 

C UVH»«VM*X{UV{3»*VMIN 

SC*L4HElGHT/aV(4>-UV(3n 

GENERAL  entry  POINT  FOR  all  AZJHOTHAl  PROJECTIONS, 

SJNPH«SIN{F*CLATI2>-PH10) » 

C6SPH«C0S(F • (LAT  r 2 ) -PH  10  > ) 

S l NRT  « -C<3S  ( F »R0  T ) 

COSRTi  3 I N ( F »RS T ) 

S l Nl0»  SINTF-FCI.AT  > 

CSSlC  *C(-)S(F»PCLAT  ) 

51NlA*SJN(F«LAT(1  ) ) 

CfcSLA»SORT(l,.SINLA»«  JNLA  ) 

CeSA«SlNLA»SINLC*COSLA*CeSL®»C0SPH 
SINa«SQRT(1,OOOO1*C0'A«COSA» 

S1NH»C0SLA»SINPH/SINA 
C6S0A(S|NLA«ceSLP-C6SLA*SlNL€*C0SPK)/3JNA 

5TERE0QRaPhIC  WITH  ORIGIN  at  PSL AT ( P0L0NQ 

Ri2,0»SINA/(1,000001«CESA) 

11(1  )■  R»(C3SB*C0SRT*SlNb»S3NRT  I 
U(2)»-R*(5lNe«C0SRT*C0SB*3INRT  I 
THIS  SECU8N  CALCULATES  SERIES  AN®  H 

I C0ULD  HAVE  BEEN  REFLACEO  PITH  Hi  BUT  IT  PAS  N0T  W0RTH  THE  EFFORT 
JjJ  IS  A COUNTER,  IT  IS  In  CEHM0N  BECAUSE  IT  K;EDS  T0  BE  INCREMENTED 
EACH  TIME  IT  SWJTCFES  FROM  «0LTPUH  T*  «C0NV» 

the  value  of  hi  is  set  in  ‘ether*  it  is  the  n^xt  unused  position 

IN  array  X <12000 ) WHICH  WILL  START  SERIES  AN8H. 

IML*K  | NE  1 1 )G6  Tfi  18 
P(1HU<1» 

P(2)«U<2> 

JJJiJJJ  ♦ 1 

K»  JjJ 
J»K-1 

If  (J)aoo, 600.31 
31  AN0M(1)«AN0M(J) 

HEDN(i )»HEDN»2) 

REA0<o6,33)ANCM(2),H6DM2HlPf’N 

33  F6RMA  T ( 2f 10 1 4 i I 2 ) 

BOO  Hill] 

IF»K,EQ,1>  501.602 
301  6 ( 1 ) »P ( 1 > 

5 < 2 J *P  < 2 > 

06  TG  20 

602  If ( I PEN , EQ , 3 > 6Q3.5Q2 

603  JPPN«3 

G6  T0  205 

502  I f ( K i EQ 1 2 ) 503,604 

503  XLIff  ■P(l)-5<n 
YClFf«P(2)«S(2) 

G6  T0  17 

604  If  UPPN,E0,3»  503,504 
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504  xcirr*(xu|FF«i»{P(i)*s<i ) ) )/(i«n 
YClFF.(rDIFF*l«(P(2>-S(2>  > > / Cl*l ) 

17  G£  TG  1? 

15  HEAD<2>»ATAn2(YD1FF,XDIFF  ) 

IF  <HEAD<2) ,CT,0, )HEAC{2)»2«PI*NfcAD(2) 

C HE  AD ( 2 ) IS  AN  ANGLE  BETWEEN  0 AND  2»FJ 

19  IF (o. , Lfc . HEAD ( 2 ) 5411,408 
All  IF (HfcAU{2) ,LE,PI/2)A07,A0B 
407  ISIgNjI 

ge  re  23 

4 08  IF  (3,  *PI/2,,LT .HEAD (2) 5412,410 
412  IF  (HfcAD{2),L7,2,*Pl>405  |410 
4Q9  IS1GN*1 
G6  7 . 23 
410  I S I QN»*1 

23  xeNE^sa  >-ISlGN»ANeM(2J*SlN(HEAD(2J)/IOINCH«SCAL) 
X7W0  *S(2)*IS1GN«ANJH<2)*CSS{H6AD(25)/IQINCH*SCAI> 
205  HH  1 tE  ( 05 , 700  5 XSNE.XTWSdPPN.lDAY 
700  F£RmAT(2F10,4,2|10) 

XL  A ST ■ XQNfc 
I HR | T« I *1 
S(1)»P(1) 

S ( 2 J »P ( 2 5 
HEAD! 1 ) (HEAD (2  5 

16  1*1*2 

20  1 I 1 1 «I 
IS  RETURN 

END 


CGNV 

KENT  C0NV 


PR0GRAP  LENGTH 

00515 

ENTRY  PeiNTS  C0NV 

00032 

6LSCK  NAMES 

00315 

1 

00004 

9 

00004 

EXTERNAL  SYMBOLS 

O1C10100 

THEND, 


Q1C04100 

Q0CCICT, 

ATAN2 

SQRTF 

SJNF 

C0SF 

TSP, 

STM , 

ONS|NOL, 

00234  SYMBOLS 
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subroutine  output  <namEiigg> 

c 

PEA).  u*T 
INTEGER  TITLE 
CIMENSleN  RAC  JUS  <2 ) 

REAL  LATMJN|LATHA*,LENGR1N,LENGPAX 
real  L*TNeT(90),LENNfT(lBO) 

CEMMGN  DEL*TiCElONiXEGLAT,XPCLGN 
C6HMGN  l» 

CGMmGn  W.DlSTjANftf-CK.CEANGEiNP.l.KK.KNUM,  INUH.QlNCH 
C6MMGN  AN0M<  2),  m.JJJ.KKK.XLASt 
COMhGN  L * T < 2 ) 

CGMmCn  JNPi1G,L,V,XLAT,XL0N 

CGMMGN  PGlATiEGL0NG,FGT,LMIN,LPAX,VmIN,VMAX,HE!GHT,NLAT,NLON, 

• IFR0J,  IBBXiSCALE.  IETPB,  JUNE 
CGMMGN  NALL,  1 start (11 >, NAMES < 10 »6>,LENG, LAST 
CGMMGN  TES»TITLEU0),lCGL»IA(7fl> 

CeMM0N/l/lDATE<2)  ,HECM2) 

CEMm0n/3/L*TPJN,LATMAX 
CEMmEn/7/LONGMN.LGNCMAX 
CGMM0N/9/UDJEF  , C IFF  ilUIN , WIN 
DATA ( F ■ ,0174533  > 

DATA <S  I N1«,017AS24)»(C€S1«J 999946) 

DATa<ENDLAT»99,0).<JFEN»0) 

C 

IHlKK.EQ,2)  GG  T0  582 

LTEmP«LKK 

LKK*C 

€NGMINsL0NGMJN 

GNGMAX«L0NGPAX 

IF(L0N«MIN,LT,O)  GNGPJN»L0NGP JN  ♦ 360 
JF<L0NGMAX,LT,O)  GNGPAX«L0NGPAX  ♦ 360 
JFEn»3 

testmjn*L0Ngm{n 
TESTMAX.L0NGMAX 
ITM]N«ABS( TESTMJN) 

ITMAX.AHS(TfcSTMAX) 

C TEST  IF  C0MPLETE  CIRCLE 

IFdTMJN. E0,0. AND, UMAX. E0J360  ) GG  J0  3000 
C TEST  IF  H0VEHS  ARGUNC  0 GR  180 

IF (TESTMlN.GT.O.ANCtTESTPAX.LT.OIGG  TB  2001 
IF (TESTMlN,LT,0,AKC,TESTPAX.GT.O)  G0  Te  2000 
IF  < ABS  <L0NGMAX  > , GT  , AES <LGNGP IN  ) )GG  T0  8500 
C LEFT  HALF  I$F  SPHERE 

IF (ABSJLQNGMlNJ.LE.OCfGR.ABSlLCNGMjNJ.OT.OO.AND.ABStLeNQMAXJ.Lt, 
190  ) GC  T 0 2005 
XL*T*LATMAX 
XL0N*  GNGMAX 
CALL  CGNV(XLAT) 
djff«v-vmin 

XL*T»LATM JN 
XL0N°0NGMAX 
CALL  CSNV(XUAT) 

UCIFFbU.UMIN 
GO  T0  8501 
2005  XLA?*l*TMJN 
XL0N*  0NGMAX 
CALL  CBNV(XLAT) 

DJFF»V'VM{N 

XL*T«L*TMIN 

XL0N»QNGMJN 

CALL  C0NVIXLAT) 

UCJFFiU.UMlN 
G0  T 0 8501 
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C COHPLETE  CJRCLE 

3000  C I F f ■ 0 
xlat*latmjn 
XL0N*270 . 0 
CALL  COnV(XlaT) 

UC  IfFxU.UMIN 

G6  TE  8531 

C ROVERS  A30UNp  180 

3001  IF  ( JTmIm.  (it,  I T MAX  ) G£  TO  2002 
XLAT«LATMAX 

XL0N»0MGMIN 

call  C0NV(XLAT) 

DIFF«V'VMJN 
XLAT»L*TMJN 
XL0NS0NGMAX 
CALL  C0NV(XLAT) 

LC  IFF  cU«UMlN 
HTEST*2 
GO  TG  05  Ji 
c hovers  around  o 
3000  DIFr»0 

XLAT»LATMJN 
XL0N»L0NGMIN 
CALL  C0NV(XLAT) 

UClrFsU-UMjN 
NLTESm 
G6  TC  85 J 1 

3003  xlat«l*tmax 

XLOfJ«0NCMAX 
CALL  C0NV(XLAT) 

D I F E «v-VMJN 
XLATtLATMIN 
XL0N«0NGMAX 
CALL  CONV(XUAT) 

UCIFFsU-UMIN 
NLTEST*2 
G6  TO  8501 

C RIGHT  HALF  OF  SPHERE 

fl9^Ol)GOAT0(2JOAMIN,,^E,95,ef,'*0S<^eNGr1*N,,l'T,  90  i *ND,  ABS(L0NGMAX) . Q T • 90 ) 00  T0 
XL*T*LATMAX 
XL0N»0NGmIN 
call  C0NV(XLAT) 

C|Ff«V*VMJN 
XLAT«LATM*X 
XL0N*0NGHAX 
call  CSNV(XLAT) 

LDJFF'U-UMIK 

G6  TO  8501 

300*  xlat«latmjn 

XL0NS0NGMJN 
CALL  C0NV(XLAT) 

CIFF«V-VHJN 

Xl*T°LATMAX 

XL0Ns0NGMJN 

CALL  CHNV(XLA1> 
lciff»u«umin 
r.e  to  05O1 
#501  VV l N* I N 
VVAx»VMAX 
UL I N»UH  I N 
ULAX*UMAX 

I C 0 U N T ■ 0 
LKK.LTEMP 
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JCHECK«0 
KEUnT.O 
LHL«1 
FIRST. 0 
100  CeNTlNUfc 
C 

C IF  •JPEf'-*  IS  0 T(- 1 S IS  THE  FIRST  HAP  AND  THE  8RIGIN  IS  N6T  SHIFTED. 

101  IF  |!PfcN>  102.103,102 

102  CALL  PLOT (WlDTH+l.CiO.C.’J) 

C 

103  WIDTH.  HE  1GHT*(IMAX-UPIM/<VPAX-VM1N) 

CALL  PLOT (0.0  .-3) 

JF(vES.tO,lO0,0)CALL  SYME0L  < *1 , 0 » . 2 . . 21 » T J TlE , 90 , 0 , 8 0 ) 

VE S. «1 1 0 

C CERTAIN  trig  FUNCTION  THAT  ARE  C6NSTANT  FOR  A GIVEN  HAP  ARE  CALCULATED  AND 
C STORED  IN  »C0NV*  * *NEWMP»  IS  AN  ENTRY  TO  THAT  ROUTINE 
C 

SCALE. HEIGHT/IVHAX* VH1N ) 

D1FF«DIFF. SCALE 
LC  IFF. UDIFF. SCALE 
LIEhPsLKK 
Lkk.q 

C CRAW  LONGITUDE  LINES 
DEG.FlOAT(NL0M 
1 F-  E n ■ 3 
VLAT'-RO.O 

PL0NG.0NGHIN 

1F(PGLONG,LT,0, )PLCNG»P0LONG436O, 

XL0N»PLONG-DEG 
DL*T  = - .5 
ZZMaX.DELAT/2,0 
ZTOp.pOLAT  ♦ ZZMAX 
DC  110  1 -1(360.  NL6N 
XL0NeXL0N*DEG 

IF(xL0N,G1 ,360, )XLeN.XL0N*36O, 
dlat«-dl*t 

DC  110  J*1 . 361 

xlat«xlat*dlat 
CALL  CONV(XLAT) 

V*MU-UMIN>*SCALE)-UCIFF 
w.( ( V-Vm I N 1 *SC ALF j-diff 
955  F6RMAT(1H0.BF1C,5) 

IF (NLTEST.EQ.I)  GO  F 767 

IF (XL0N.LT,  ONGMIN..*,  L€N:GT,  eNGHAX,0R.XLAT,LT.LATMIN,0H,XLAT, 
1GT.LATHAX)  GC  TO  120 
GC  TO  778 

767  lF(XL0N,GT,0NGHAX.ANC,XLeNlLT,GNGHlN,8R,XLAT.LT,LATMlH,OR,XLAT,OT, 
1LATMAX)  GO  TO  120 
776  IF ( XL AT. ZT0P., A >7000.120.7000 
7000  CALL  plot  ( V . W i I PEN  ) 

IPEN*2 
G6  T0  110 
120  I P E N ■ 3 

110  continue 

DEG.Fl0aT(NlAT  I 
150  XLAT.LATMIN-CEG 
IC2.2.NLAT 

c 

ce  lt>9  I * I G2  , J61 , NL  AT 
xl*t«xlat*deg 

IF (XLAT.GT ,L*THAX  ♦ 1)  G«  Tl  666 
IF (XLAT , GE , 90 , >900. SOI 
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900  XLAT»L*TM jn.CEG 
CE0»"DtG 

901  DL0N*1 . C 
I FEn*  3 
Xw0N"-Dl0N 

155  XLRN«XL0N«DL6N 
CALL  CBNV(XLAT) 

W«<  (V.VMlN)*StAl,E)-DlFF 
VM(ll-UMlN>»SrALF!-LClFF 
IF (NGTEST,EU,1)  G£  TG  779 
556  IF  (XC0N.LT  , GNGM  I N , (S  P , XLGN1  GT  , 
1 G 7 .LATMAX)  OG  T0  156 
GE  TC  780 


6hCPAX,0R.XUT,LT.lATMlN,0H,X|.AT, 


779 

780 

156 

158 

159 
C 


1LATHA®?'?0  '?SG^8  * ANC  ,XUeN*LT  *0NGMIN',R*  KLAT  'LT*LA™IN*C,,,Xl*AT»aT* 

CALL  PL0T(V,W, IFEN) 

IPEn  * 2 

IF { XL  0M- 360, 0)155, 159, 159 
IPEn* J 
G£  TO  156 
CENT INUE 


C LABEL  LATITUDE  LINES 
886  I5T0P.LONQMAX  . LEKGPJN  ♦ 1 
I T Op  ■ UEl*T  *2*1 
IPEn  « J 
DEGsF  lOaT (NlAT ) 

IF (NUTEbT ,EQ,1)  G£  T £ 775 
XPOLON*  OnGMA*.(  £M,MX»  8KGMNI/2 
G£  TO  774 

775  XPO[>0N»L0NGHAX*(L£NGPAX*LeisGMM/2 
774  xpolat«latmin  - DEG 
D6  171  I *1 • 1 TCP 
XP0LAT»XPBLAT  * DEG 
IF(XP0LAT,GT,LATMAX)  GO  T0  999 
CALL  C^NV ( XPELAT ) 

V*((L.UMIN)*SCALE)-LTIFF 
w*( ( V* Vm I N ) • S C A L 6 ) -D I FF 
CALL  PL0T(Y.,C5,W« ,05,3> 

171  CeNTINUE,9fctt<V4,05't,4'05,  ,07,)<FetA1'  0t0*4HF6,l) 

c label  longitude  lines 

999  DEGcFLOAT  (NLGM 

XP0LAT«L*TMJN  « DELAT/2 
XP0LON»L0NGMIN  - DEG 
609  XP0lONsXP0LON  * DEG 
A * 27 0 , ♦ XP0LEN 
IF(a,GT,360J  A«a-360 
IF  < NL1TE5T  , Eu, 2 ) Gg  1 £ 6503 
IF (xP0LON,GT,L0NG6AX)  GO  TE  998 
GE  TO  85u A 

8503  IF(XP0L0N,KT,GNPVAX)  GO  Tg  998 
650«  CALL  CONV(XPGLAT) 

W*( (V.VMIN)»SC*LE)-CIFF 
V*( (1-UMIN)*SCALE).IEJFF 
CALL  PLOT ( 7*, Q5 , w* , 05  »3  ) 

CAL^.  NOMBEP(Y*,05tl»-,05i  , 07 1 XP6L6N , A» 4HF6 , 1 ) 

Cg  TO  8o5 
996  LKK'LTEMP 

C 

c 
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C 

REWJND  20 
582  11*1 

IEND»l.*ST  *2 

IF i m» j end  • 4 

J y U o 0 

|FEN*3 

INMAPrO 

c the  value  of  lkk  is  set  in  slb  ether  depending  bn  the  value  of  kkk 

IF  < LkK  . fc  Q , 2 ) GE  TO  J03 
GE  T 0 301 

309  H»<  (LAM2)-VHIM«SCALb)-C!FF 
Y*<  <LAT<l)-UMIN)*SCALfc)-lD!FF 
GE  TE  310 

C 

C CHECK  F0R  BEGINNING  £F  NEW  SERIES  6R  CHANGE  T0  POINT  M0DE 
301  REAd(20,6000)LAT<1)iLAT(2),IPEN  ,IPATE(2> 

6000  F£RMAT<2F10, 4,2110) 

IF (LAT<l)-ENCLAT)3C2,328,32e 
C 

C CONTINUOUS  M0DE  CAT*  DRAWN 
302  CALL  C0NV(LAT(  1 ) i I PEN , 1 PR£y ) 

GE  TO  304 

303  RE  AD ( 05 , 6000  ' LA  ‘ ( 1 > , LA  T < 2 ) , I PEN  » 1 DATE  < 2 > 

I F (EOF i 05)390*7) 

777  IF  <LAT<1>»ENCIAT)309.329i30« 

329  INMaP*INMAP*1 
IPEn-3 
n*i  i*2 

331  IF  ( I I-IMN)303,390.350 
304  CONTINUE 

W«( (V.VMIN)*SCALE)-DIFF 
Y»< (U*UMIN)*SCALE)»UCIFF 
C 

C CHECKS  IP  POINT  LIES  INS  ICE  MAP  RECTANGLE.  IF  NOT  SKIPS  PLOT  ROUTINE  AND 

c counting  statement 

310  I F ( I PEN , EQ , 5 ) GO  T6  210 
9067  CALL  PL0T(Y,W,IPEN) 

210  CENTlNUfc 

INMaP«InMaP*1 
I PEn*2 
320  11-11*2 

IF<LKK.EQ,2)G6  TO  331 
Ic  ( I J - I 6ND  ) 301. 390. 390 
328  INMAP«lNMAP*l 
330  IFEn»3 

GE  TO  320 

390  I H *P * JNMAp  - l 

IF (LKK.EQ.2)  CO  T£  604 
WRITE( 10.3000)  IMAP 

3000  permatiiohoprcghah  plotted. 19, 2xi7Hp0|nts  on  the  map) 

C IF  LKK«0  SU‘J  OTHER  wll  NET  CALCULATE  ANOMALY  SERIES  { IF»8  ANOMALY  SERIES 
C MAC  JUST  BEEN  PLOTTED  THEREFOR  REINITIALIZE  EVERYTHIN!) 

IF (lKK.EO.O.OR.LKK.EC ,2)604 »399 
6Q4  KNUm»0 
INUm«0 
Uv  JaO 
NALL*0 

last«o 

399  CONTINUE 


28 


o o c* 


NRL  REPORT  7929 


c pregham  plotted  track  ready  ti  Pier 

REWJND  05 
500  RETURN 

PRINTED  OUTPUT  SECTION 
END 


ICENT 


PREGRAM  LENGTH 

0 2161 

ENTRY  P6INTS 

OUTPUT 

0 0 < 5 6 

BLOCK  NAMES 

1 

0C315 
0 0 0 0 4 

3 

00002 

7 

00002 

9 

00004 

external  symbols 

01C10100 

THEND, 

QHGE  ICT. 

CONV 

PLET 

SYPEEL 

NUMBER 

O8CIPE0F 

REW, 

TSP, 

STH, 

ONSiNGL, 

00423  SYMBOLS 


ANOMALY 


output 
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SCBReuTjNfc  NEXT(lNSTR,NAKE> 

C 

c this  program  attempts  te  frcvice  a machine-independent  rsutine  res  reading 

C CQNTRLu  CARDS  IN  SUPERMAP,  THE  L£RC  LENGTH  0F  THE  MACRJNE  IS  REQUIRED  T0  BE 
C AT  LEAST  FuUR  BCD  CHARACTERS  L£NG.  THIS  IS  MET  BY  ALL  MACHINES  LIKELY  V0  BE 
C ENCELNTERED. 

c 

DIMENSION  N AM’E  ( 6 ) 

CEMMON  DELAT,CEL0N,XF0LAT,XPSLPN 
CEMmCn  lx 

CGMM0N  W.DIST, ANOMCK I CHANGE . N P . L*K , KNUM , INUM.GINCH 
CEMmEn  AN0M  ( 2),IIJ,Jjj,Kt<K,XLAST 
C&MmSN  X ( 2 ) 

CEMmFN  I NP  i 1 0 , L) , V , XL  * T i XL0N 

CEMHEN  M0L*T iFELONG (R8T,LMIN,LMAX,VmJN,VMAX,HEIGHT,NLAT|NL0Ni 
♦IFR0J, IB0X(SCAL6. ISYMB,ILINE 
CEMHEN  NALL, ! START (11) .NAMES UO|0> i LENG»LAST 
CEMHEN  YES  » T I TLE  ( 10  ) I JCGL  # I A < 76 ) 

cemmgn/i/ i date (2 t » he dn ( 2 ) 

data ( j blank oih  > . < jcgmm»ih,  > , (Mll«1) 

10  CE  11  I W 0 R D * 1 , 6 

11  NAME ( I WORD ) ■ IELANk 
J KQrD*1 

15  IF  ( I W0RD i GT 1 1 ) RETURN 
20  JC0L«JCOL*1 

IF  ( I COL  < GT , 76  ) G0  TS  3Q 
25  J A C « I A ( I CQL ) 

IF  ( I AC , E 0 . IBL ANK  ,OR,  J AC jEC  . ICGMMa)  G0  T9  15 
IF  ( IW0HD < LE , 6 ) NAME( JL0RD)«IA< IC6L) 

IKORC«IWORn*l 

null»i 

GE  TE  20 

30  IF  (NULL  tEQ, 0 ) G0  TE  35 

IF  ( I HERD , G T , 1 1 RETURN 
READ  (IMP,  1000  ) I NS  TP , I A 
1000  FERMAT(AA,76A1) 

NULL'O 

NAME(1)«IBLANK 
D6  33  KcOL«1i76 

C KC6L  IS  USED  BECAUSE  THE  C0MPILEF!  SEEMS  UNABLE  T0  ACCEPT  THE  09  LMP  BEL9H  IF 
C I C6L  IS  USED  THREUQHeiT, 

IC0L»KC0L 

IF  ( lA(KC8L),EQ,IBLANK.eR,  IA(KC6L)  ,EO. JC0MMA)  Q0  T0  10 
33  CENT INUE 

35  NULL*! 

RETURN 

END 


PR8GRAM  LENGTH 
ENTRY  PUNTS  NEXT 
0L9CK  NAMES 

1 

EXTERNAL  SYMB0LS 

ThEND, 

Q9CCICT, 

TSH, 

SLl. 

QN3JNGL, 

00190  5TMB0LS 


IDEM 

00200 

00012 

00315 

00004 


NEXT 
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900 

200 

eoe 


201 

977 

eoo 


779 

777 

776 


SLBHBUT  j NE  QTKRJXELCjXNEW) 

RfcAt  LATHlN,L*TMAX,LEKGMlN,i.eNOPAX 

120 >^AAioMAU  JsSI20  * *,,CAY<20,*t  <20  > ,DDMIN<20  > . R R l.  A T ( 2 0 > . RRL 0N Q ( 

DIMENSION  IFM  (20  > •KRFEM2  ) 

CEHMPN  UfcLAT  .CELO\,  XFQ(.A7  , yFELGN 
CEMmKn  LX 

C6HH0N  ^nMfT,AN^C^^A,NGEuK1,L,K'KNUM*INUH'GlNCH 
CfcMHON  AN0M(  2 ) I I II  I Jv J, KKK  , XLASt 

CEMH0N  X ( 2 ) 

CEHMEN  I NP , IE,U,V,XLAT,XL0K 

CtMH0N/l/IDATE(2)  ,HEDM2) 

CEHM0N/3/LATMIN,LATPAX 
CEMH0N/5/JUDY1, JUCY2, JTMl, MP2 
C£Mm0m/7/l0nGPIN«LCNGMAX 

CEHM0M/e/jrM 

CEMmGn/U/IeXTRAi  ISKJP 

ITSJ5 
I E X * 0 

atmjn.latmin 

AYMAX.LATMAX 

I F ( J S YHfcj . fcO . 1 ) A Tp  J N» ATM  J N ♦ 1 
D l S T 2 « G 0 . 

HEDrl.279, 

I FL  J GHT»1 
I K P « 6 0 

IAY1«3 

icay»o 

ZN I NE  «99 . 0 
ZE«0«O  . 0 

CEGrA.1,745J2SE-2 
CMNGE1.45°,,(3ftO,-ChAAGE/2, ) 

CEAnGE2B(450(0 -CHANGE/2 , )-J80, 

P * 1 
L»l 

IF  <LKK.EQ,1>G£  TO  50 
I«1 

IF  <l.E3, 1 )200, 201 
WMTE(2C.701)ZMNE,ZERE,IDAY,ICAY 
FEMmATCIHO*  HO,  2F10.4) 

P«M*2 

N*2l 

IF  (N-20)  801,601,977 
Nil 

XA*NgmIl[n^M>I,YR<N>,IICAY<N>'1 IHR<K  > »DdH*NCN> ,RRLAt(N) ,RRU9N0(N* , 

I F t I ECMfcCK , I T ) 800,779 
IF(Eer.jT)  777,778 
IENDMEND  * 1 

I F < J END  . GE,  1)  G6  T8  401 
N»N  • 1 

IP<N*lMl>  06  TO  600 
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N»1 

001  I Y G * I 1 Y H ( \ ) 

!LAr*IlDAY<N» 

1 ♦-  R a I JHh(N) 

CF I Nl  * bDM  1 Nl  ( N ) 

JE  INSOMJN 

!KIN«1HW*100  * JMIN 

RL*T'HRl*T ( N ) 

RL0NG«RNl_0Nli<N> 

*N0MAL*AAN0MAL (N  ) 

IF ( IB©X 13^1,355,350 
391  AN0MAL=»*N0MAb 
GG  T 0 355 

390  an0maub»an0mac»i . e 2 e e 

399  N«N*1 

IF ( 1 n A Y ,11 , JUL v 1 ) G€  1G  201 
IF (RLAT  ,GT  ,V0.0)O£  10  202 
IF ( I C A Y , G T , JILY2)  00  TC  Ell 
IF ( I Y R . 1 0 . 0 ) 00  H Bll 
IF ( RL A T ,LT  , ATKIN)  G0  T0  202 
IF  <RUT,GT,LATPAX>  <>E  10  202 
IF (RU0NU.UT,UeKGMlN)  GO  10  202 
IF (HL0NG.GT  .LGKGMAX ) G0  10  202 
IF  ( JDAY  ,EQ,  JUCY1.  AND,  IE  JN.U  • HE1  > GO  T0  201 
IF < jDaY.EQ, JUCY2. AND. IE JN.G1 , I1E2)  GO  T0  811 
IF ( IER  .fcQ  , EHR , AND, WE IN  ,U1 .KE IN ) GG  JO  202 
401  IF <c-3>70|72.72 
202  I * Y l * 3 

GG  T0  201 

811  IF<f3>V>,  810.610 

810  Rfi  I TE  < 06 . 71 ) *NGK  < 1 ) . F-EE N < 1 ) » HPP6N  < 1 J 
GE  T0  9V 

72  RR  I TE  ( 06  • 71 ) ANGM(1),HEDM1)iHPPEM1) 

71  FERMAT! 2F 10, 4|I2) 

70  IF ( lfcND,EO,l>  GE  IE  09 
501  IF  <L,fcQ,l)GO  K 5 
RED'H  1 ) *MfcCM 
9 GGNTlNUfc 

AN0M(1)«AN0M(2) 

KFPEN(l)*KPPENf2) 

KFP£N ( 2 ) * I a Y1 

AN0M(2)«AN0MAL 

GG  T0  19 
19  cgntinue 

RF.lTE(20.70l)fiLAT,RL€NGiIAYl»IAYl 

701  FGRMATI2F10, 4.2J10) 

KF  R * I HR 
KF IN’JMJN 
I A Y 1 s 2 

IF ( t , EQ , 1 ) G0  T6  60 
MEDn<  I)p450,-HECM  II 

IF(HEDM(  1),GT,30O.OO>HEDM  1)pHEDN(  D-360.0 
910  HEDN(  lltHEDM  1 ) *CEGRA 
60  CGNTINUE 
l«l*l 

H«M*2 
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I f < J t G T , 1 0 > UO  T0  2C 
12  KK  = J *1 
?0  I'M 

ce  T6  50 J 

99  WRlTE(06>71)AN£H(2)iHEEN(l)|KPPt;Nt2) 

I E X. I fcX  ♦ 1 

irtjExTKA.EO, JEX)  G6  TS  667 
L«1 

P»|T  ♦ 1 

READ(6C,665>ISKJPiwUCV1, J TW1 , UICY2 , I TH2 
669  F E KMA  T (9l4) 

DE  781  I X«l( ISK 1 P 
781  CALL  SKJPMtECIT) 

J A Y1 *3 
I END  ■ 0 
G£  TG  977 
667  LENg*m-1 
RbwjND  06 
I I I «H 

M*l-1 

IF (KKK.EQ,0)206i207 
206  LKK.Q 
RETURN 


207  LKK.l 
RETURN 

rorTrE  V*LUE  0E  HI  IE  SET  IN  5L8  CGNV  WHEN  SERIES  TRACK  IS  BEING  PL0TTED 
0*  EAch  !f*LUfc  0F  TRACK  PL0TTEL  A CGRRE  SPGND l NQ  VALUE  0F  SERIES  AN6M 

(The  series  will  el  an  x»y  value  in  inches  stsredin  common  xciodoi 
khen  series  an9e  is  flatted  SLBfieLUNE  cbnvert  is  bypassed  X(ICOO) 

SINCE  THE  SERIES  AN0F  IS  ALREADY  IN  INCHES 
III  IS  0NE  MSRE  TE*N  TEE  TElAL  LENTH  0F  aHRAYX  ( 10000  ) 

LENGIh  m each  series  calculated  in  *0ther«  which  is  stbred 
lenqI^i ^1  ^3  **N6fIj**  THE  LENGTH  G,r  THE  used  OBRTIBN  er  array  xuoooo) 

LKKS2  * * 

RETgRN 
END 


90 


PR8GRAH  LENGTH 

ICENT 

01049 

entry  peints 

0THER 

0C242 

BLICK  NAMES 

1 

00315 
0000  4 

3 

00002 

9 

0CCC4 

7 

.00002 

8 

00024 

10 

00002 

EXTERNAL  SYMB0LS 

00330  SYMB0LS 
•binary  deck**# 

THEND, 

OICIOIOO 

Q6CUCT, 

SKIFEILE 

Q8CIFE0F 

Q8CJFI0C 

REW, 

TSH, 

onsIngl. 

BANK, <0)|/1/ 
L0AD 

RUN, 90, 10000 


•ther 
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